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We present ten new equilibrium solutions to plane Couette flow in small periodic cells 
at low Reynolds number Re and two new traveling- wave solutions. The solutions are 
continued under changes of Re and spanwise period. We provide a partial classification 
of the isotropy groups of plane Couette flow and show which kinds of solutions are 
allowed by each isotropy group. We find two complementary visualizations particularly 
revealing. Suitably chosen sections of their 3D-physical space velocity fields are helpful in 
developing physical intuition about coherent structures observed in low Re turbulence. 
Projections of these solutions and their unstable manifolds from their oo-dimensional 
state space onto suitably chosen 2- or 3-dimensional subspaces reveal their interrelations 
and the role they play in organizing turbulence in wall-bounded shear flows. 

1. Introduction 

In Gibson et al. (2008) (henceforth referred to as ) we formed visualizations of 
the oo-dimensional state-space dynamics of moderate Re turbulent flows, using precisely 
calculated equilibrium solutions of the Navier- Stokes equations to define dynamically in- 
variant, intrinsic, and representation independent coordinate frames. These state-space 
portraits (figure 1) offer a visualization of the dynamics of transitionally turbulent flows, 
complementary to 3D visualizations of the spatial features of velocity fields (figure 2). 
Side-by-side animations of the two visualizations illustrate their complementary strengths 
(see Gibson (20086) online simulations). In these animations, 3D spatial visualization of 
instantaneous velocity fields helps elucidate the physical processes underlying the forma- 
tion of unstable coherent structures, such as the Self-Sustained Process (SSP) theory of 
Waleffe (1990, 1995, 1997). Running concurrently, the oo-dimensional state- space repre- 
sentation enables us to track the unstable manifolds of equilibria and the heteroclinic 
connections between them (Halcrow et al. 2009), and provides us with new insight into 
the nonlinear state space geometry and dynamics of moderate Re wall-bounded flows. 

Here we continue our investigation of equilibrium and traveling-wave solutions of 
Navier-Stokes equations, presenting ten new equilibrium solutions and two new traveling 
waves of plane Couette flow, and continuing the solutions as functions of Re and periodic 
cell size [L X ,2,L Z ]. Nagata found the first pair of nontrivial equilibria (Nagata 1990) 
and the first traveling wave in plane Couette flow (Nagata 1997). Clever & Busse (1992) 
found closely related equilibria in plane Couette flow with Rayleigh-Benard convection. 
Cherhabili & Ehrenstein (1997) reported two-dimensional equilibria of plane Couette 
but these were later shown to be artifacts of the truncation (Rincon 2007; Ehrenstein 
et al 2008). Waleffe (1998, 2003) computed the Nagata equilibria guided by the SSP 
theory and showed that they were insensitive to the boundary conditions at the wall. 
Other traveling waves were computed by Viswanath (2008) and Jimenez et al (2005). 
Schmiegel (1999) computed and investigated a large number of equilibria. His 1999 Ph.D. 
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Figure 1. A 3- dimensional projection of the oo-dimensional state space of plane Couette flow 
in the periodic cell Q G hc at Re = 400, showing all equilibria and traveling waves discussed in 
§4. Equilibria are marked: EQ (laminar flow), O EQ 1? • EQ 2 , □ EQ 3 , ■ EQ 4 , EQ 5 , < 
EQ 7 , < EQ 8 , EQ 9 , V EQ 10 , T EQ n , A EQ 12 , and A EQ 13 . Traveling waves trace out closed 
orbits: the spanwise-traveling TWi (blue loops), streamwise TW2 (green lines), and TW3 (red 
lines). In this projection the latter two streamwise traveling waves appear as line segments. The 
EQi — > EQ heteroclinic connections and the S- invariant portion of the EQ X and EQ 2 unstable 
manifolds are shown with black lines. The cloud of dots are temporally equispaced points on a 
long transiently turbulent trajectory, indicating the natural measure. The projection is onto the 
translational basis (3.13) constructed from equilibrium EQ 2 . 

provides a wealth of ideas and information on solutions to plane Couette flow, and in 
many regards the published literature is still catching up this work. added the dy- 

namically important unb equilibrium (labeled EQ 4 in this paper). Parallel theoretical 
advances have been made in channel and pipe flows, with the discovery of traveling waves 
in channel flows (Waleffe 2001), and traveling waves (Faisst & Eckhardt 2003; Wedin & 
Kerswell 2004; Pringle & Kerswell 2007) and relative periodic orbits (Duguet et at. 2008) 
in pipes. Moreover, traveling waves have been observed experimentally in turbulent pipe 
flow (Hof et al. 2004). We refer the reader to GHC for a more detailed review. 

We review plane Couette flow in § 2. The main advances reported in this paper are 
(§ 3) a classification of plane Couette symmetry groups that support equilibria, (§ 4) the 
determination of a number of new equilibria and traveling waves of plane Couette flow, 
and (§5, §6) continuation of these solutions in Reynolds number and spanwise aspect 
ratio. Outstanding challenges are discussed in § 7. Detailed numerical results such as 
stability eigenvalues and symmetries of corresponding eigenfunctions are given in Halcrow 
(2008), while the complete data sets for the invariant solutions can be downloaded from 
channelf low . org. 

2. Plane Couette flow — a review 

Plane Couette flow is comprised of an incompressible viscous fluid confined between 
two infinite parallel plates moving in opposite directions at constant velocities, with no- 
slip boundary conditions imposed at the walls. The plates move along the streamwise 
or x direction, the wall-normal direction is y, and the spanwise direction is z. The fluid 
velocity field is u(x) = [u,v,w](x,y, z). We define the Reynolds number as Re = Uh/v, 
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Figure 2. A snapshot of a typical turbulent state in a large aspect-ratio cell 
[L X ,2,L Z ] = [15,2,15], Re = 400. The walls at y = ±1 move away/towards the viewer at 
equal and opposite velocities U = ±1. The color indicates the streamwise (u, or x direction) 
velocity of the fluid: red shows fluid moving at u — +1, blue, at u — —1. The colormap as 
a function of u is indicated by the laminar equilibrium in figure 4. Arrows indicate in-plane 
velocity in the respective planes: [v, w] in (y, z) planes, etc. The top half of the fluid is cut away 
to show the [u,w] velocity in the y = midplane. See Gibson (20086) for movies of the time 
evolution of such states. 

where U is half the relative velocity of the plates, h is half the distance between the plates, 
and v is the kinematic viscosity. After non-dimensionalization, the plates are positioned 
at y = ±1 and move with velocities u = ±lx, and the Navier-Stokes equations are 

?^ + u • Vu = - Vp + -J- V 2 u , V • u = . (2.1) 
at Re 

We seek spatially periodic equilibrium and traveling- wave solutions to (2.1) for the do- 
main Q = [0, L x ] x [—1, 1] x [0, L z ] (or Q = [L x , 2, L z ]), with periodic boundary conditions 
in x and z. Equivalently, the spatial periodicity of solutions can be specified in terms of 
their fundamental wavenumbers a and 7. A given solution is compatible with a given 
domain if a = m2n/ L x and 7 = n2n/L z for integer m, n. In this study the spatial mean 
of the pressure gradient is held fixed at zero. 

Most of this study is conducted at Re — 400 in one of the two small aspect-ratio cells 

Q Gnc = [2tt/1.14, 2, 2tt/2.5] w [5.51, 2, 2.51] w [190, 68, 86] wall units 

Q HKW = [2tt/1.14, 2, 2tt/1.67] w [5.51, 2, 3.76] w [190, 68, 128] wall units, (2.2) 

where the wall units are in relation to a mean shear rate of (du/dy) = 2.9 in non- 
dimensionalized units computed for a large aspect-ratio simulation at Re = 400. Empir- 
ically, at this Reynolds number the HKW cell sustains turbulence for very long times 
(Hamilton et al. 1995), whereas the GHC cell exhibits only short-lived transient tur- 
bulence (GHC). The z length scale L z = 4tt/5 of fl GHC was chosen as a compromise 
between the L z = 6tt/5 of H kw 

and its first harmonic L z /2 = 3tt/5 (Waleffe 2002). 
Unless stated otherwise, all calculations are carried out for Re = 400 and the ^ GHC cell. 
In the notation of this paper, the solutions presented in Nagata (1990) have wavenum- 
bers (a, 7) = (0.8,1.5) and fit in the cell [2tt/0.8, 2, 2tt/1.5] w [7.85, 2, 4.18]. f Waleffe 
(2003) showed that these solutions first appear at critical Reynolds number of 127.7 and 

f Note also that Reynolds number in Nagata (1990) is based on the full wall separation and 
the relative wall velocity, making it a factor of four larger than the Reynolds number used in 
this paper. 
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(a, 7) = (0.577,1.15). Schmiegel (1999)'s study of plane Couette solutions and their 
bifurcations was conducted in the cell of size = [471, 2, 2n] « [12.57, 2, 6.28]. 

Although the aspect ratios studied in this paper are small, the 3D states explored 
by equilibria and their unstable manifolds explored here are strikingly similar to typi- 
cal states in larger aspect-ratio cells, such as figure 2. Kim et al (1971) observed that 
streamwise instabilities give rise to pairwise counter-rotating rolls whose spanwise sepa- 
ration is approximately 100 wall units. These rolls, in turn, generate streamwise streaks 
of high and low speed fluid, by convecting fluid alternately away from and towards the 
walls. The streaks have streamwise instabilities whose length scale is roughly twice the 
roll separation. These 'coherent structures' are prominent in numerical and experimental 
observations (see figure 2 and Gibson (20086) animations), and they motivate our inves- 
tigation of how equilibrium and traveling-wave solutions of Navier-Stokes change with 
Re and cell size. 

Fluid states are characterized by their energy E = ^||u|| 2 and energy dissipation rate 
D = || V x u|| 2 , defined in terms of the inner product and norm 

(u,v) = i /dxu-v, ||u|| 2 = (u,u). (2.3) 

The rate of energy input is I = 1/(L X L Z ) J J dxdzdu/dy, where the integral is taken 
over the upper and lower walls at y = ±1. Normalization of these quantities is set so that 
/ = D = 1 for laminar flow and E = / — D. It is often convenient to consider fields as 
differences from the laminar flow, since these differences constitute a vector space, and 
thus can be added together, multiplied by scalars, etc. We indicate such differences with 
tildes: u = u — yx.. Note that the total velocity field u does not form a vector space: 
the sum of any two total plane Couette velocity fields violates the u = ±1 boundary 
conditions at the moving walls. 



3. Symmetries and isotropy subgroups 

On an infinite domain and in the absence of boundary conditions, the Navier-Stokes 
equations are equivariant under any 3D translation, 3D rotation, and x — > — x, u — > — u 
inversion through the origin (Frisch 1996). In plane Couette flow, the counter- moving 
walls restrict the rotation symmetry to rotation by tt about the z-axis. We denote this 
rotation by cr x and the inversion through the origin by <j xz . The suffixes indicate which 
of the homogeneous directions x, z change sign and simplify the notation for the group 
algebra of rotation, inversion, and translations presented in §3.1 and §3.2. The cr xz and 
a x symmetries generate a discrete dihedral group D\ x D\ = {e, a x , cr z , a xz } of order 4, 
where 

a x [u, v, w] (£, y, z) = [-u, -v, w](-x, -y, z) 

a z \u,v,w\(x,y,z) = [u,v,-w](x,y,-z) (3.1) 
a xz [u,v,w](x,y,z) = [-u, -v, -w](-x, -y, -z) . 

The walls also restrict the translation symmetry to 2D in-plane translations. With 
periodic boundary conditions, these translations become the SO (2) x x SO(2) z continuous 
two-parameter group of streamwise-spanwise translations 

r(£ x ,£ z )[u,v,w](x,y,z) = [u, v, w]{x + 4, y, z + t z ) . (3.2) 

The equations of plane Couette flow are thus equivariant under the group T = 0(2) x x 
0(2) z = Di x x SO(2) x x D\ z x S0(2) Z1 where x stands for a semi-direct product, x 
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subscripts indicate streamwise translations and sign changes in x,y, and z subscripts 
indicate spanwise translations and sign changes in z. 

The solutions of an equivariant system can satisfy all of the system's symmetries, a 
proper subgroup of them, or have no symmetry at all. For a given solution u, the subgroup 
that contains all symmetries that fix u (that satisfy su = u) is called the isotropy (or 
stabilizer) subgroup of u. (Hoyle 2006; Marsden & Ratiu 1999; Golubitsky & Stewart 
2002; Gilmore & Letellier 2007). For example, a typical turbulent trajectory u(x, t) has 
no symmetry beyond the identity, so its isotropy group is {e}. At the other extreme is 
the laminar equilibrium, whose isotropy group is the full plane Couette symmetry group 

r. 

In between, the isotropy subgroup of the Nagata equilibria and most of the equilibria 
reported here is S = {e, si, 82, S3}, where 

si [u,v,w](x,y,z) = [u,v,-w](x + L x /2,y,-z) 

s 2 [u, v, w](x, y, z) = [-u, -v, w](-x + L x /2, -y, z + L z /2) (3.3) 
s 3 [u, v, w](x, y, z) = [-u, -v, -w](-x, -y, -z + L z /2) . 

These particular combinations of flips and shifts match the symmetries of instabilities of 
streamwise-constant streaky flow (Waleffe 1997, 2003) and are well suited to the wavy 
streamwise streaks observable in figure 2, with suitable choice of L x and L z . But S is 
one choice among a number of intermediate isotropy groups of T, and other subgroups 
might also play an important role in the turbulent dynamics. In this section we provide a 
partial classification of the isotropy groups of T, sufficient to classify all currently known 
invariant solutions and to guide the search for new solutions with other symmetries. We 
focus on isotropy groups involving at most half-cell shifts. The main result is that among 
these, up to conjugacy in spatial translation, there are only five isotropy groups in which 
we should expect to find equilibria. 

3.1. Flips and half-shifts 

A few observations will be useful in what follows. First, we note the key role played by 
the rotation and reflection symmetries cr x and <j z (3.1) in the classification of solutions 
and their isotropy groups. The equivariance of plane Couette flow under continuous 
translations allows for traveling- wave solutions, i.e., solutions that are steady in a frame 
moving with a constant velocity in (x,z). In state space, traveling waves either trace 
out circles or wind around tori, and these sets are both continuous-translation and time 
invariant. The sign changes under <j x , <j z , and <j xz , however, imply particular centers of 
symmetry in x, z : and both x and z, respectively, and thus fix the translational phases 
of fields that are fixed by these symmetries. Thus the presence of a x or a z in an isotropy 
group prohibits traveling waves in x or z, and the presence of a xz prohibits any form 
of traveling wave. Guided by this observation, we will seek equilibria only for isotropy 
subgroups that contain the cr xz inversion symmetry. 

Second, the periodic boundary conditions impose discrete translation symmetries of 
r(L x ,0) and r(0,L z ) on velocity fields. In addition to this full-period translation sym- 
metry, a solution can also be fixed under a rational translation, such as r(mL x /n, 0) or a 
continuous translation r(4,0). If a field is fixed under continuous translation, it is con- 
stant along the given spatial variable. If it is fixed under rational translation r(mL x /n^ 0), 
it is fixed under r(mL x /n, 0) for m G [1, n — 1] as well, provided that m and n are rela- 
tively prime. For this reason the subgroups of the continuous translation SO(2) x consist 
of the discrete cyclic groups C UyX for n = 2,3,4,... together with the trivial subgroup 
{e} and the full group SO(2) x itself, and similarly for z. For rational shifts £ X /L X = m/n 
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we simplify the notation a bit by rewriting (3.2) as 

r^ n = r(mL x /n, 0) , r^ n = r(0, mL z /n) . (3.4) 
Since m/n = 1/2 will loom large in what follows, we omit exponents of 1/2: 

If a field u is fixed under a rational shift r(L x /n), it is periodic on the smaller spatial 
domain x G [Q,L x /n]. For this reason we can exclude from our searches all equilibrium 
whose isotropy subgroups contain rational translations in favor of equilibria computed on 
smaller domains. However, as we need to study bifurcations into states with wavelengths 
longer than the initial state, the linear stability computations need to be carried out 
in the full [L X ,2,L Z ] cell. For example, if EQ is an equilibrium solution in the = 
[L X /3,2,L Z ] cell, we refer to the same solution repeated thrice in Q = [L X ,2,L Z ] as 
"spanwise-tripled" or 3 x EQ. Such solution is by construction fixed under the Cs yX = 
{e,ri /3 ,rj /3 } subgroup. 

Third, some isotropy groups are conjugate to each other under symmetries of the full 
group T. Subgroup H' is conjugate to H if there is an s G T for which H' = s~ 1 Hs. 
In spatial terms, two conjugate isotropy groups are equivalent to each other under a 
coordinate transformation. A set of conjugate isotropy groups forms a conjugacy class. 
It is necessary to consider only a single representative of each conjugacy class; solutions 
belonging to conjugate isotropy groups can be generated by applying the symmetry 
operation of the conjugacy. 

In the present case conjugacies under spatial translation symmetries are particularly 
important. Note that 0(2) is not an abelian group, since reflections a and translations 
r along the same axis do not commute (Harter 1993). Instead we have err = r~ l cr. 
Rewriting this relation as err 2 = r _1 <jr, we note that 

(4,0) = r" 1 (4/2,0) 

0~x T\t x AO). (3.6) 

The right-hand side of (3.6) is a similarity transformation that translates the origin of 
coordinate system. For £ x = L x /2 we have 

T- 1/4 ^ry 4 = <r x r x , (3.7) 

and similarly for the spanwise shifts / reflections. Thus for each isotropy group containing 
the shift-reflect cr x r x symmetry, there is a simpler conjugate isotropy group in which cr x r x 
is replaced by <j x (and similarly for cr z r z and a z ). We choose as the representative of each 
conjugacy class the simplest isotropy group, in which all such reductions have been made. 
However, if an isotropy group contains both cr x and cr x r xi it cannot be simplified this 
way, since the conjugacy simply interchanges the elements. 

Fourth, for t x = L x , we have r" 1 cr x r x = cr x , so that, in the special case of half-cell 
shifts, cr x and r x commute. For the same reason, cr z and r z commute, so the order- 16 
isotropy subgroup 

G = D liX x C 2 , x x D ltZ x C 2 , z C T (3.8) 

is abelian. 

3.2. The 67-fold path 

We now undertake a partial classification of the lattice of isotropy subgroups of plane 
Couette flow. We focus on isotropy groups involving at most half-cell shifts, with SO(2) x x 
SO(2) z translations restricted to order 4 subgroup of spanwise-streamwise translations 



Equilibria and traveling waves of plane Couette flow 



7 



(3.5) of half the cell length, 

T = C 2 , x x C 2iZ = {e,r x ,T z ,T xz } . (3.9) 

All such isotropy subgroups of V are contained in the subgroup G (3.8). Within G, we 
look for the simplest representative of each conjugacy class, as described above. 

Let us first enumerate all subgroups H C G. The subgroups can be of order \H\ = 
{1,2,4,8, 16}. A subgroup is generated by multiplication of a set of generator elements, 
with the choice of generator elements unique up to a permutation of subgroup elements. 
A subgroup of order \H\ = 2 has only one generator, since every group element is its own 
inverse. There are 15 non-identity elements in G to choose from, so there are 15 subgroups 
of order 2. Subgroups of order 4 are generated by multiplication of two group elements. 
There are 15 choices for the first and 14 choices for the second. However, each order-4 
subgroup can be generated by 3 • 2 different choices of generators. For example, any two 
of t x1 t z1 t xz in any order generate the same group T. Thus there are (15 • 14) /(3 • 2) = 35 
subgroups of order 4. 

Subgroups of order 8 have three generators. There are 15 choices for the first generator, 
14 for the second, and 12 for the third. There are 12 choices for the third generator and 
not 13, since if it were chosen to be the product of the first two generators, we would 
get a subgroup of order 4. Each order-8 subgroup can be generated by 7 • 6 • 4 different 
choices of three generators, so there are (15 • 14 • 12)/(7 • 6 • 4) = 15 subgroups of order 
8. In summary: there is the group G itself, of order 16, 15 subgroups of order 8, 35 of 
order 4, 15 of order 2, and 1 (the identity) of order 1, or 67 subgroups in all (Halcrow 
2008). This is whole lot of isotropy subgroups to juggle; fortunately, the observations of 
§3.1 show that there are only 5 distinct conjugacy classes in which we can expect to find 
equilibria. 

The 15 order-2 groups fall into 8 distinct conjugacy classes, under conjugacies between 
a x r x and a x and a z r z and a z . These conjugacy classes are represented by the 8 isotropy 
groups generated individually by the 8 generators a x , o Zl o~ xz , o~ x t Zi o~ z t Xi r X: r z , and 
r xz . Of these, the latter three imply periodicity on smaller domains. Of the remaining 
five, a x and cr x r z allow traveling waves in z, a z and cr z r x allow traveling waves in x. Only 
a single conjugacy class, represented by the isotropy group 

{e,(T xz }, (3.10) 

breaks both continuous translation symmetries. Thus, of all order-2 isotropy groups, we 
expect only this group to have equilibria. EQ 9 , EQ 10 , and EQ n described below are 
examples of equilibria with isotropy group {e,a xz }. 

Of the 35 subgroups of order 4, we need to identify those that contain <j xz and thus 
support equilibria. We choose as the simplest representative of each conjugacy class the 
isotropy group in which a xz appears in isolation. Four isotropy subgroups of order 4 
are generated by picking a xz as the first generator, and o Zl o~ z t Xi o~ z t Zi or g z t xz as the 
second generator (R for reflect-rotate) : 

R = {e, a x , a z , a xz } = {e, a xz } x {e, a z } 

R x = {e,o- x T x ,o- z T x ,a xz } = {e,a xz } x {e,a x r x } (3.11) 
R z = {e, a x T z , cr z T z , a xz } = {e, a xz } x {e, a z r z } 

Rxz ~ Ox^xz-) ^z^xz-i ®~xz\ ~ ®~xz\ X {fi, O ' Z T XZ \ — S . 

These are the only isotropy groups of order 4 containing a xz and no isolated translation 
elements. Together with {e, cr xz }, these 5 isotropy subgroups represent the 5 conjugacy 
classes in which expect to find equilibria. 
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The R xz isotropy subgroup is particularly important, as the Nagata (1990) equilibria 
belong to this conjugacy class (Waleffe 1997; Clever & Busse 1997; Waleffe 2003), as do 
most of the solutions reported here. The NBC isotropy subgroup of Schmiegel (1999) and 
S of Gibson et al. (2008) are conjugate to R xz under quarter-cell coordinate transforma- 
tions. In keeping with previous literature, we often represent this conjugacy class with 
S {e,si,s 2 ,s 3 } = {e, 

o~z7~xi o~xi~xzi o~xzi~z\ rather than the simpler conjugate group R xz . 
Schmiegel's I isotropy group is conjugate to our R z \ Schmiegel (1999) contains many 
examples of R z -isotropic equilibria. i?-isotropic equilibria were found by Tuckerman & 
Barkley (2002) for plane Couette flow in which the translation symmetries were broken 
by a streamwise ribbon. We have not searched for ^-isotropic solutions, and are not 
aware of any published in the literature. 

The remaining subgroups of orders 4 and 8 all involve {e,Ti} factors and thus involve 
states that are periodic on half-domains. For example, the isotropy subgroup of EQ 7 and 
EQ 8 studied below is S x {e,r xz } ~ R x {e,r xz }, and thus these are doubled states of 
solutions on half-domains. For the detailed count of all 67 subgroups, see Halcrow (2008). 



3.3. State-space visualization 

GHC presents a method for visualizing low-dimensional projections of trajectories in the 
infinite-dimensional state space of the Navier-Stokes equations. Briefly, we construct an 
orthonormal basis {ei,e 2 , • • • ,e n } that spans a set of physically important fluid states 
ua,U5,..., such as equilibrium states and their eigenvectors, and we project the evolving 
fluid state u(t) = u(t) — yx onto this basis using the L 2 inner product (2.3). This produces 
a low-dimensional projection 

a(t) = (ai, a 2 , • • • , a n , • • • )(t) , a n (t) = (u(t), e n ) , (3.12) 

which can be viewed in 2d planes {e m ,e n } or in 3d perspective views {e^, 

^n}« The 

state-space portraits are dynamically intrinsic, since the projections are defined in terms 
of intrinsic solutions of the equations of motion, and representation independent, since 
the inner product (2.3) projection is independent of the numerical or experimental repre- 
sentation of the fluid state data. Such bases are effective because moderate- Re turbulence 
explores a small repertoire of unstable coherent structures (rolls, streaks, their mergers), 
so that the trajectory a(t) does not stray far from the subspace spanned by the key 
structures. 

There is no a priori prescription for picking a 'good' set of basis fluid states, and 
construction of {e n } set requires some experimentation. Let the ^-invariant subspace be 
the flow- invariant subspace of states u that are fixed under S; this consists of all states 
whose isotropy group is S or contains S as a subgroup. The plane Couette system at 
hand has a total of 29 known equilibria within the /S-invariant subspace: four translated 
copies each of EQ X - EQ 6 , two translated copies of EQ 7 (which have an additional r xz 
symmetry), plus the laminar equilibrium EQ at the origin. As shown in GHC, the 
dynamics of different regions of state space can be elucidated by projections onto basis 
sets constructed from combinations of equilibria and their eigenvectors. 

In this paper we present global views of all invariant solutions in terms of the or- 
thonormal 'translational basis' constructed in GHC from the four translated copies of 
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Figure 3. Four projections of equilibria, traveling waves and their half-cell shifts onto transla- 
tional basis (3.13) constructed from equilibrium EQ 4 . Equilibria are marked EQ , O EQ l5 • 
EQ 2 , □ EQ 3 , ■ EQ 4 , EQ 5 , < EQ 7 , EQ 9 , V EQ 10 , T EQ n , A EQ 12 , and A EQ 13 . Traveling 
waves trace out closed loops. In some projections the loops appear as line segments or points. 
TWi(blue) is a spanwise- traveling, symmetry-breaking bifurcation off EQ 1? so it passes close to 
different translational phases of O EQ X . Similarly, TW3 (red) bifurcates off □ EQ 3 and so passes 
near its translations. TW2 (green) was not discovered through bifurcation (see §4); it appears 
as the shorter, isolated line segment in (ai,a4) and (02,04). The EQ X — > EQ relaminarizing 
heteroclinic connections are marked by dashed lines. A long-lived transiently turbulent trajec- 
tory is plotted with a dotted line. The EQ 4 -translational basis was chosen here since it displays 
the shape of traveling waves more clearly than the projection on the EQ 2 -translational basis of 
figure 1. 
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(3.13) 



where c n is a normalization constant determined by ||e n || = 1. The last 3 columns indicate 
the symmetry of the basis vector under half-cell translations; e.g. ±1 in the r x column 
implies r x ej = ±ej. 



4. Equilibria and traveling waves of plane Couette flow 

We seek equilibrium solutions to (2.1) of the form u(x, t) = u EQ (x) and traveling- wave 
or relative equilibrium solutions of the form u(x,t) = u TW (x — ct) with c = (c x ,0,c z ). 
Let Fns( u ) represent the Navier-Stokes equations (2.1) for the given geometry, boundary 
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conditions, and Reynolds number, and f^ s its time-t forward map 

^=F NS (u), f^ s (u) = u + jf drP NS (u(r)). (4.1) 

Then for any fixed T > 0, equilibria satisfy f T (u) — u = and traveling waves satisfy 
f T (u) — r u = 0, where r = r(c x T, c z T). When u is approximated with a finite spectral 
expansion and f £ with a numerical simulation algorithm, these equations become a set 
of nonlinear equations in the expansion coefficients for u and, in the case of traveling 
waves, the wave velocities (c x ,0,c z ). 

Viswanath (2007) presents an algorithm for computing solutions to these equations 
based on Newton search, Krylov subspace methods, and an adaptive 'hookstep' trust- 
region limitation to the Newton steps. This algorithm can provide highly accurate so- 
lutions from even poor initial guesses. The high accuracy stems from the use of Krylov 
subspace methods, which can be efficient with 10 5 or more spectral expansion coeffi- 
cients. The robustness with respect to initial guess stems from the hookstep algorithm. 
The hookstep limitation restricts steps to a radius r of estimated validity for the local lin- 
ear approximation to the Newton equations. As r increases from zero, the direction of the 
hookstep varies smoothly from the gradient of the residual within the Krylov subspace 
to the Newton step, so that the hookstep algorithm behaves as a gradient descent when 
far away from a solution and as the Newton method when near, thus greatly increasing 
the algorithm's region of convergence around solutions, compared to the Newton method 
(J.E. Dennis, Jr., & Schnabel 1996). 

The choice of initial guesses for the search algorithm is one of the main differences 
between this study and previous calculations of equilibria and traveling waves of shear 
flows. Prior studies have used homotopy, that is, starting from a solution to a closely 
related problem and following it through small steps in parameter space to the problem of 
interest. Equilibria for plane Couette flow have been continued from Taylor- Couette flow 
(Nagata 1990), Rayleigh-Benard flow (Clever & Busse 1997), and from plane Couette with 
imposed body forces (Waleffe 1998). Equilibria and traveling waves have also been found 
using "edge-tracking" algorithms, that is, by adjusting the magnitude of a perturbation 
of the laminar flow until it neither decays to laminar nor grows to turbulence, but instead 
converges toward a nearby weakly unstable solution (Itano & Toh (2001); Skufca et al. 
(2006); Viswanath (2008); Schneider et al. (2008)). In this study, we take as initial guesses 
samples of velocity fields generated by long-time simulations of turbulent dynamics. The 
intent is to find the dynamically most important solutions, by sampling the turbulent 
flow's natural measure. 

We discretize u with a spectral expansion of the form 

J K L 

u ( x )=E E J2 u ^ T ^ e2Mix/L * +kz/Lz) > ( 4 - 2 ) 

j = -J k=-K £=0 

where the are Chebyshev polynomials. Time integration of f £ is performed with a 
primitive-variables Chebyshev-tau algorithm with tau correction, influence-matrix en- 
forcement of boundary conditions, and third-order semi-implicit backwards-diflerentiation 
time stepping, and dealiasing in x and z (Kleiser & Schumann (1980); Canuto et al. 
(1988); Peyret (2002)). We eliminate from the search space the linearly dependent spec- 
tral coefficients of u that arise from incompressibility, boundary conditions, and com- 
plex conjugacies that arise from the real-valuedness of velocity fields. Our codes for 
Navier-Stokes integration, Newton-hookstep search, parametric continuation, and eigen- 
value calculation are available for download from channelflow.org website (Gibson 
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2008a), along with a database of all solutions described in this paper. For further details 
on the numerical methods see GHC and Halcrow (2008). 

Solutions presented in this paper use spatial discretization (4.2) with (J, K, L) = 
(15, 15,32) (or 32x33x32 gridpoints) and roughly 60k expansion coefficients, and inte- 
gration is performed with a time step of At = 0.03125. The estimated accuracy of each 
solution is listed in table 1. As is clear from Schmiegel (1999) Ph.D. thesis, ours is al- 
most certainly an incomplete inventory; while for any finite i?e, finite aspect-ratio cell the 
number of distinct equilibrium and traveling wave solutions is finite, we know of no way 
of determining or bounding this number. It is difficult to compare our solutions directly 
to those of Schmiegel since those solutions were computed in a [4tt, 2, 2n] cell (roughly 
twice our cell size in both span and streamwise directions) and with lower spatial reso- 
lution (2212 independent expansion functions versus our 60k for a cell of one- fourth the 
volume) . 

4.1. Equilibrium solutions 

Our primary focus is on the ^-invariant subspace (3.3) of the ^ GHC cell at Re = 400. We 
initiated 28 equilibrium searches at evenly spaced intervals At = 25 along a trajectory in 
the unstable manifold of EQ 4 that exhibited turbulent dynamics for 800 nondimension- 
alized time units after leaving the neighborhood of EQ 4 and before decaying to laminar 
flow. Lower, upper branch pairs are labeled with consecutive numbers, and the numbers 
indicate, as closely as possible, the order of discovery. We give a name EQ n to any distinct 
solution in the ^ GH c cell at Re = 400, although many of these solutions can be connected 
by continuation in Re or wavenumber. EQ is the laminar equilibrium, EQ X and EQ 2 
are the Nagata lower and upper branch, and EQ 4 is the u NB solution reported 
The rest are new. Only one of the 28 searches failed to converge onto an equilibrium; the 
successful searches converged to equilibria with frequencies listed in table 1. The higher 
frequency of occurrence of EQ X and EQ 4 suggests that these are the dynamically most 
important equilibria in the ^-invariant subspace for the ft GHC cell at Re = 400. Stability 
eigenvalues of known equilibria are plotted in figure 7. Tables of stability eigenvalues and 
other properties of these solutions are given in Halcrow (2008), while the images, movies 
and full data sets are available online at channelflow.org. All equilibrium solutions 
have zero spatial-mean pressure gradient, which was imposed in the flow conditions, and, 
due to their symmetry, zero mean velocity. 

EQ l5 EQ 2 equilibria. This pair of solutions was discovered by Nagata (1990), re- 
computed by different methods by Clever & Busse (1992, 1997) and Waleffe (1998, 
2003), and found multiple times in searches initiated from turbulent simulation data, as 
described above. The lower branch EQ X and the upper branch EQ 2 are born together 
in a saddle-node bifurcation at Re ~ 218.5. Just above bifurcation, the two equilibria 
are connected by a EQ X — > EQ 2 heteroclinic connection, see Halcrow et al. (2009). How- 
ever, at higher values of Re there appears to be no such simple connection. The lower 
branch EQ 1 equilibrium is discussed in detail in Wang et al (2007). This equilibrium has 
a 1-dimensional unstable manifold for a wide range of parameters. Its stable manifold 
appears to provide a partial barrier between the basin of attraction of the laminar state 
and turbulent states (Schneider et al. 2008). The upper branch EQ 2 has an 8-dimensional 
unstable manifold and a dissipation rate that is higher than the turbulent mean, see fig- 
ure 8(a). However, within the S- invariant subspace EQ 2 has just one pair of unstable 
complex eigenvalues. The two-dimensional /S-invariant section of its unstable manifold 
was explored in some detail in GHC. It appears to bracket the upper end of turbulence 
in state space, as illustrated by figure 1. 
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Figure 4. Equilibrium solutions of plane Couette flow in Q G hc = [27t/1.14, 2, 2tt/2.5] at 
Re = 400. Plotting conventions are the same as figure 2. The color scale indicates the streamwise 
velocity u, with the front face of the laminar solution EQ , u(y) = y, serving as a reference. Not 
all solutions are shown; EQ 3 , EQ 6 , EQ 7 , EQ 12 , and EQ 13 , are very similar to EQ 4 , EQ 5 , EQ 8 , 
EQ 10 , and EQ n , respectively. 

EQ 3 , EQ 4 . The upper branch EQ 4 solution was found in GHC and is called u NB there. 
Its lower-branch partner EQ 3 was found by continuing EQ 4 downwards in Re and also 
by independent searches from samples of turbulent data. EQ 4 is, with EQ X , the most 
frequently found equilibrium, which attests to its importance in turbulent dynamics. 
Like EQ l5 EQ 4 serves as a gatekeeper between turbulent flow and the laminar basin 
of attraction. As shown in GHC, there is a heteroclinic connection from EQ 4 to EQ X 
resulting from a complex instability of EQ 4 . Trajectories on one side of the heteroclinic 
connection decay rapidly to laminar flow; those one the other side take excursion towards 
turbulence. 

EQ 5 , EQ 6 . The lower-branch EQ 5 solution was found only once in our 28 searches, 
and its upper-branch partner EQ 6 only by continuation in Reynolds number. We were 
only able to continue EQ 6 up to Re = 335. At this value it is highly unstable, with a 
19 dimensional unstable manifold, and it is far more dissipative than a typical turbulent 
trajectory. 

EQ 7 , EQ 8 appear together in a saddle node bifurcation in Re (see §5). EQ 7 / EQ 8 
might be the same as Schmiegel (1999) 's L a solutions'. The x-average velocity field plots 
appear very similar, as do the D versus Re bifurcation diagrams. We were not able to 
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Figure 5. x-average of difference from laminar flow, u = u — yx, for equilibria EC^-EQg in 
the Qghc cell. The axes are z (horizontal) and y (vertical). Arrows indicate [v,^], with the 
same arrow length denoting the same magnitude of in-plane velocity in all graphs. The colormap 
indicates u: red/blue is u = ±1, and green is u = 0. Lower and upper-branch pairs are grouped 
together vertically, e.g. EQ l5 EQ 2 are a lower, upper branch pair. Re = 400 except for Re = 330 
in EQ 6 . 




Figure 6. x-average of u for equilibria EQ 7 -EQ 11 in the Q G hc cell at Re = 400. Plotting 
conventions are the same as in figure 5. 



obtain Schmiegel's data in order to make a direct comparison. EQ 7 is both the closest 
state to laminar in terms of disturbance energy and the lowest in terms of drag. It has 
one strongly unstable real eigenvalue within the ^-invariant subspace and two weakly 
unstable eigenvalues with {51,53} and {52,53} antisymmetries, respectively. In this re- 
gard, the EQ7 unstable manifold might, like the unstable manifold of EQ l5 form part 
of the boundary between the laminar basin of attraction and turbulence. EQ 7 and EQ 8 
are unique among the equilibria determined here in that they have the order-8 isotropy 
subgroup S x {e, r xz } (see § 3.2). The action of the quotient group Gj (S x {e, r xz }) yields 
2 copies of each, plotted in figure 3. EQ 7 and EQ 8 are similar in appearance to EQ 5 and 
EQ 6 , except for the additional symmetry. 

EQ 9 is a single lopsided roll-streak pair. It is produced by a pitchfork bifurcation 
from EQ 4 at Re « 370 as an {si, ^-antisymmetric eigenfunction goes through marginal 
stability (the only pitchfork bifurcation we have yet found) and remains close to EQ 4 
at Re = 400. Thus, it has {e,<j xz } isotropy. Even though EQ 9 is not ^-isotropic, we 
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Figure 7. Eigenvalues of equilibria EQ 3 , EQ 4 and EQ 7 , EQ 8 in the Q G hc cell, Re 



400. 



Eigenvalues are plotted according to their symmetries: • + + +, the ^-invariant subspace, ► 

H , M — I — , and A h, where zb symbols stand for symmetric/antisymmetric in si, S2, 

and S3 respectively. For EQ l5 EQ 2 and EQ 4 eigenvalues see (there referred to as u LB , u UB 

and u N b, respectively). For numerical values of all stability eigenvalues see Halcrow (2008) and 
channelf low . org. 



(a) 
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(b) 
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Figure 8. Rate of energy input at the walls / versus dissipation D, for equilibria in the (a) 
1]ghc (b) Qhkw cell at Re = 400. Typical turbulent trajectories marked are with dotted lines. 
EQ is the laminar equilibrium with D = I = 1. In (a), O EQ l5 □ EQ 3 , ■ EQ 4 , < EQ 7 , and 
EQ 9 , are clustered in the range 1.25 < D — I < 1.55, and • EQ 2 , EQ 5 , A EQ 8 , V EQ 10 , T 
EQ n A EQ 12 and A EQ 13 lie in 2 < D = I < 4. Traveling waves TWi, TW 2 , and TW 3 (not 
shown) are clustered together with O EQ l5 etc. In (b), the symbols are the same, with both the 
lower and upper branches of EQ 9 marked with . ■ EQ 4 , • EQ 2 are clustered together near 
D = I w 2.47. 
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Table 1. Properties of equilibrium solutions for Qghc cell, Re — 400, unless noted otherwise. The 
mean values are ensemble and time averages over transient turbulence. ||-|| is the L 2 -norm of the 
velocity deviation from laminar, E is the energy density (2.3), D is the dissipation rate, H is the 
isotropy subgroup, dim W u is the dimension of the equilibrium's unstable manifold or the number 
of its unstable eigenvalues, and dimVK^ is the dimensionality of the unstable manifold within 
the i^-invariant subspace, or the number of unstable eigenvalues with the same symmetries 
as the equilibrium. The accuracy acc. of the solution at a given resolution (a 32 x 33 x 32 
grid) is estimated by the magnitude of the residual ||(/ T=1 (u) — u)||/||u|| when the solution is 
interpolated and integrated at higher resolution (a 48x49x48 grid, At = 0.02). The freq. column 
shows how many times a solution was found among the 28 searches initiated with samples of 
the natural measure within the ^-invariant subspace. See also figure 8 (a). 



found it from a search initiated on a guess that was ^-isotropic to single precision. Such 
small asymmetries were enough to draw the Newton-hookstep search algorithm out of 
the ^-invariant subspace. 

EQ 10 , EQ n are produced in a saddle- node bifurcation at Re « 348 as a lower / up- 
per branch pair, and they lie close to the center of mass of the turbulent repeller, see 
figure 8(a). The velocity fields have a similar appearance to typical turbulent states 
for this cell size. However, they are both highly unstable and unlikely to be revisited 
frequently by a generic turbulent fluid state. Their isotropy subgroup {e,a xz } is order 
2, so the action of the quotient group G/{e, r xz } yields 8 copies of each, which appear 
as the 4 overlaid pairs in projections onto the (3.13) basis set, see figure 1 and figure 3. 
EQ 12 , EQ 13 were found by continuation of EQ 10 in 7. 

4.2. Traveling waves 

The first two traveling-wave solutions reported in the literature were found by Nagata 
(1997) by continuing EQ X equilibrium to a combined Couette / Poiseuille channel flow, 
and then continuing back to plane Couette flow. The result was a pair of streamwise trav- 
eling waves arising from a saddle-node bifurcation. Viswanath (2008) found two traveling 
waves, 'Dl' and the same solution 'D2,' but at a higher Re = 1000, through an edge- 
tracking algorithm (Skufca et al (2006), see also §4.1). Here we verify Viswanath's 
solution and present two new traveling-wave solutions computed as symmetry-breaking 
bifurcations off equilibrium solutions. We were not able to compare these to Nagata's 
traveling waves since the data is not available. The traveling waves are shown as 3D 
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Figure 9. Spanwise TWi, streamwise TW2 and TW3 traveling waves in Q G hc cell, Re = 400. 



|| -|| E D H dim W u dim W% c mean u 



mean 0.2828 0.087 2.926 

TWi 0.2214 0.1341 1.510 {e,a x r z } 3 2 0.00655 z 0.00482 z 

TW 2 0.1776 0.1533 1.306 {e,a z r x } 3 2 0.3959 x 0.0879 x 

TW 3 0.2515 0.1520 1.534 {e,a z r x } 4 2 0.4646 x 0.1532 x 

Table 2. Properties of traveling- wave solutions for (] G hc cell, Re — 400, denned as in table 1, 
with wave velocity c and mean velocity. See also figure 8 (a). 



velocity fields in figure 9 and as closed orbits in state space in figure 3. Their kinetic 
energies and dissipation rates are tabulated in table 1. Each traveling- wave solution has 
a zero spatial-mean pressure gradient but non-zero mean velocity in the same direction 
as the wave velocity. It is likely each solution could be continued to zero wave velocity 
but non-zero spatial- mean pressure gradient. 

TWi traveling wave is ^-isotropic and hence spanwise traveling. At Re = 400 its 
velocity is very small, c = 0.00655 z, and it has a small but nonzero mean velocity, 
also in the spanwise direction. This is a curious property: TWi induces bulk transport 
of fluid without a pressure gradient, and in a direction orthogonal to the motion of the 
walls. TWi was found as a pitchfork bifurcation from EQ l5 and thus lies very close to it 
in state space. It is weakly unstable, with a 3D unstable manifold with two eigenvalues 
which are extremely close to marginal. In this sense TWi unstable manifold is nearly 
one-dimensional, and comparable to EQ X . 

TW2 is a streamwise traveling wave found by Viswanath (2008) and called Dl there. 
It is si -isotropic, has a low dissipation rate, and a small but nonzero mean velocity in the 
streamwise direction. Viswanath provided data for this solution; we verified it with an 
independent numerical integrator and continued the solution to ^ GH c cell f° r comparison 
with the other traveling waves. In this cell TW2 is fairly stable, with an eigenspectrum 
similar to TWi's, except with different symmetries. 

TW3 is an si -isotropic streamwise traveling wave with a relatively high wave velocity 
c = 0.465 x and a nonzero mean velocity in the streamwise direction. Its dissipation rate 
and energy norm are close to those of TWi • 

5. Continuation under Reynolds number 

The relations between the equilibrium and traveling-wave solutions can be clarified by 
tracking their properties under changes in Re and spatial periodicity. Figure 10 shows 
a bifurcation diagram for equilibria and traveling waves in the GHC cell, with dissipa- 
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(a) R e (b) Re 

Figure 10. (a) Dissipation D of equilibria as a function of Reynolds number for the (] G hc cell, 
a, 7 = 1.14,2.5. O RQ,, • EQ 2 , □ EQ 3 , ■ EQ 4 , EQ 5 , EQ 6 , < EQ 7 , < EQ 8 , * EQ 9 , V 
EQ 10 , T EQ n , AEQ 12 , AEQ 13 , and > TWi. (b) Detail of (a). EQ l5 EQ 2 , EQ 3 , EQ 4 , EQ 7 , and 
EQ 9 extend past Re = 1000 (not plotted here). 



tion rate D plotted against Re as the bifurcation parameter. A number of independent 
solution curves are shown in superposition. This is a 2-dimensional projection from the 
oo-dimensional state space, thus, unless noted otherwise, the apparent intersections of 
the solution curves do not represent bifurcations; rather, each curve is a family of so- 
lutions with an upper and lower branch, beginning with a saddle-node bifurcation at a 
critical Reynolds number. 

The first saddle-node bifurcation gives birth to the Nagata lower branch EQ X and upper 
branch EQ 2 equilibria, at Re ~ 218.5. EQ 1 has a single S'-isotropic unstable eigenvalue 
(and additional subharmonic instabilities that break the S'-isotropy of EQ X ). Shortly after 
bifurcation, EQ 2 has an unstable complex eigenvalue pair within the ^-invariant subspace 
and two unstable real eigenvalues leading out of that space. As indicated by the gentle 
slopes of their bifurcation curves, the Nagata (1990) solutions are robust with respect to 
Reynolds number. The lower branch solution has been continued past Re = 10, 000 and 
has a single unstable eigenvalue throughout this range (Wang et al. 2007). 

EQ 7 and EQ 8 are formed in a saddle-node bifurcation at Re ~ 235. The EQ 8 bi- 
furcation curve is unusual in that it increases rapidly from the bifurcation point to a 
maximum dissipation of D = 5.2 at Re = 364, and then turns rapidly but smoothly 
back down to much lower dissipation at higher Reynolds numbers. This behavior persists 
when examined at higher spatial and temporal resolutions. 

EQ 3 and EQ 4 were discovered in independent Newton searches and subsequently found 
by continuation to be lower and upper branches of a saddle-node bifurcation occurring 
at Re « 364. EQ 3 has a leading unstable complex eigenvalue pair within the S'-invariant 
subspace. Its remaining two unstable eigen-directions are nearly marginal and lead out 
of this space. 

EQ 6 was found by continuing EQ 5 backwards in Re around the bifurcation point at 
Re « 326. We were not able to continue EQ 6 past Re = 335. At this point it has a 
nearly marginal stable pair of eigenvectors whose isotropy group is T, which rules out a 
bifurcation to traveling waves along these modes. Just beyond Re = 335 the dynamics in 
this region appears to be roughly periodic, suggesting that EQ 6 undergoes a supercritical 
Hopf bifurcation here. At Re ~ 348, EQ 10 , EQ n are born in a saddle node bifurcation, 
similar in character to the EQ X / EQ 2 bifurcation, as are EQ 12 , EQ 13 . 

Figure 10(b) shows several symmetry-breaking bifurcations. At Re ~ 250, TWi bi- 
furcates from EQ X in a subcritical pitchfork as an ^-symmetric, si, ^-antisymmetric 
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eigenfunction of EQ X becomes unstable, resulting in a spanwise- moving traveling wave. 
At Re ~ 370, the EQ 9 equilibrium bifurcates off EQ 4 along an si, s 2 -antisymmetric, 
53-symmetric eigenfunction of EQ 4 . Since S3 symmetry fixes phase in both x and z, this 
solution bifurcates off EQ 4 as an equilibrium rather than a traveling wave. 

Figure 11(a) shows the instability of each equilibrium as a function of Reynolds num- 
ber. As a measure of instability we use the sum of the real parts of the equilibrium's 
unstable eigenvalues, i.e. the local exponential rate of stretching of the equilibrium's 
unstable manifold. Several observations can be made. Each lower-branch solution (open 
symbol) is less unstable than its upper-branch counterpart (closed symbol) . Lower branch 
solutions tend towards lesser instability as Reynolds number increases, upper branch so- 
lutions become more unstable. Since lower/upper branches of a given solution are defined 
by lower/higher dissipation rates, this implies that lower instability and lower dissipation 
go hand in hand. However this relation does not generally hold between different solution 
branches: EQ 7 has lower dissipation than EQ X (figure 10(a)) but is more unstable (Fig- 
ure 11(a)). Several of the lower-branch solutions have very slowly decreasing instability 
over the range of Reynolds numbers shown; for these, the number of unstable eigenvalues 
is constant or slowly decreasing as well. EQ X has three unstable eigenvalues shortly after 
bifurcation and just one for 270 < Re < 10, 000, EQ 7 has six after bifurcation and three 
for 340 < Re < 800, EQ 3 has four from bifurcation onwards, 363.9 < Re < 800, and EQ 12 
has eight from bifurcation onwards, 324.4 < Re < 600. The upper limits of these ranges 
are merely the endpoints of our calculations. Clever & Busse (1997) show that at certain 
wavenumbers, EQ 2 is stable for a small range of Reynolds numbers just after bifurcation. 
We did not look for or find regions of stability for any of the new solutions, though we 
expect such regions exist for carefully tuned parameters. For example, a stable periodic 
orbit has been found in the related Kuramoto-Sivashinsky system (Lan & Cvitanovic 
2008). For the upper branch solutions generally both the numbers of unstable eigenval- 
ues as well as their sums increase with Re. Lastly, it should be remembered that it is 
not at all clear how much the instability of an equilibrium has to do with the Lyapunov 
exponents of the turbulent flow - that depends on how close and how frequently a typical 
trajectory visits the neighborhood of a given equilibrium. 

6. Continuation under spanwise wavenumber 

In this section we examine changes in solutions under variation in spanwise period- 
icity. Figure 11(b) shows the dissipation D of the solutions as a function of spanwise 
wavenumber 7. Only a few of the intersections in this plot indicate bifurcations; the rest 
are artifacts of the projection onto the (j : D) plane. The true bifurcations are EQ n 
branching off from EQ 2 near (7, D) = (2, 3); EQ 3 from EQ 7 near (7, D) = (2.2, 1.2); and 
EQ 9 from EQ 4 near (7, D) = (2.6, 1.4). Continuation in 7 also shows that the EQ 10 , EQ n 
and EQ 12 , EQ 13 solution curves, which appear to be independent in figure 10, are con- 
nected by a saddle-node bifurcation between EQ 10 and EQ 12 near (7, £>) = (2.15,2.05). 

In addition to these bifurcations, we are interested in connecting the solutions for 
^ghc discussed in § 4, to the wider ft HKW cell of Hamilton et al. (1995), which empirically 
exhibits turbulence for long time scales at Re = 400. Of the equilibria discussed above, 
only EQ 4 , EQ 7 , EQ 8 , and EQ 9 could be continued at Re = 400 from 7 = 2.5 of H GHC 
down to the fundamental wavenumber 7 = 1.67 of ^ H kw EQ 9 appears in a saddle- 
node bifurcation just below 7 = 1.67. The other equilibria terminate in saddle- node 
bifurcations above 7 = 1.67 or at bifurcations from other solution curves. However, the 
EQ X and EQ 2 solutions can be continued upwards in 7 to the first harmonic 7 = 3.34 = 
2 x tt/L z of ft 

hkw (L z — 1.67). These solutions then appear in Q HKW as as spanwise 
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Figure 11. (a) Sum of the real parts of unstable eigenvalues of equilibria as a function of 
Re, at a,7 = 1.14,2.5 (the Q GH c cell), (b) Dissipation D of equilibria as a function of spanwise 
wavenumber 7, with fixed streamwise wavenumber a = 1.14 and Reynolds number Re = 400. 
Symbols are O EQ l5 • EQ 2 , □ EQ 3 , ■ EQ 4 , EQ 5 , ♦ EQ 6 , < EQ 7 < EQ 8 , EQ 9 , V 
EQ 10 , T EQ n . A EQ 12 , A EQ 13 . The vertical dotted line marks the fundamental wavenumber 
7 = 1.67 = 2n/L z of the Q H kw cell (L z — 3.76), the dashed line the fundamental 7 = 2.5 of Q G hc 
(L z — 2.51), and the dot-dashed line the first harmonic 7 = 3.34 of Q H kw- The intersection of 
the EQ l5 EQ 2 curves with 7 = 3.34 indicates that these solutions exist as "doubled" states in the 
Dhkw cell, as shown in figure 12. The EQ 7 and EQ 8 curves join smoothly near (7, D) — (2.6, 5.4). 
We were not able to continue EQ 6 past 7 = 2.38, nor ■ EQ 4 below 7 = 1.64. 




"doubled" states 2 x EQ X and 2 x EQ 2 . Figure 12 shows 3D velocity fields for equilibria 
in D H k W : EQ 4 , EQ 7 , EQ 8 , and EQ 9 and the spanwise doubled 2 x EQ X and 2 x EQ 2 . 
The upper branch of EQ 9 not shown, is very similar to EQ 9 . The properties of these 
solutions are listed in table 6. 

The low dissipation values of the HKW equilibria in figure 8 (b) suggest that they 
are not involved in turbulent dynamics, except perhaps as gatekeepers to the laminar 
equilibrium. We suspect that the equilibria as yet undiscovered, or the already known 
periodic orbit solutions (Gibson & Cvitanovic 2009) do play a key role in organizing 
turbulent dynamics. However, unlike the GHC cell, we were not able to find any equilibria 
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Table 3. Properties of equilibrium and traveling wave solutions for Q H kw cell, Re = 400, 
denned as in table 1. For 2xEQ x and 2xEQ 2 , H is the 8th-order group generated by t z ,ct z Tx, 
and GxztV^ . See also figure 8 (b). 



for [] H kw cell from initial guesses sampled from turbulent trajectories within the S- 
invariant subspace. This is curious, contrasted to our success in finding equilibria from 
such guesses in £l GH c? and it suggests that the aspect ratios of the HKW cell are the 
most incommensurate (fit the intrinsic widths of rolls least well) compared to the roll 
and streak scales of spanwise-infinite domains, which are apparent (approximately) in the 
simulation of figure 2. The stability calculations by Clever & Busse (1997) indicate that 
the Nagata solutions prefer a 2:1 streamwise to spanwise aspect ratio. Hence a study of 
changes in solutions under variation in both streamwise and spanwise periodicities might 
shed further light on the physical nature of these solutions. 

7. Conclusion and perspectives 

As a turbulent flow evolves, every so often we catch a glimpse of a familiar pattern. For 
any finite spatial resolution, the flow approximately follows for a finite time a pattern 
belonging to a finite alphabet of admissible fluid states, represented here by a set of 
equilibrium and traveling wave solutions of Navier-Stokes. These are not the 'modes' of 
the fluid; they do not provide a decomposition of the flow into a sum of components at 
different wavelengths, or a projection basis for low-dimensional modeling. Each solution 
spans the whole range of physical scales of the turbulent fluid, from the outer wall-to-wall 
scale, down to the viscous dissipation scale. Numerical computations require sufficient 
resolution to cover all of these scales, so no global dimension reduction is likely. The role 
of invariant solutions of Navier-Stokes is, instead, to partition the oo-dimensional state 
space into a finite set of neighborhoods visited by a typical long-time turbulent fluid 
state. 

Motivated by the recent observations of recurrent coherent states in experiments and 
numerical studies, we undertook here an exploration of the hierarchy of all known equilib- 
ria and traveling waves of fully-resolved plane Couette flow in order to describe the spatio- 
temporally chaotic dynamics of transitionally turbulent fluid flows. Turbulent plane Cou- 
ette dynamics visualized in state space appears pieced together from close visitations to 
coherent states connected by transient interludes, as can be seen in Gibson (20086) ani- 
mations of figure 2. The 3D fluid states explored by the small aspect-ratio equilibria and 
their unstable manifolds studied in this paper are strikingly similar to states observed in 
larger aspect-ratio simulations, such as figure 2. 
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For plane Couette flow equilibria, traveling waves and periodic solutions embody a 
vision of turbulence as a repertoire of recurrent spatio-temporal patterns explored by 
turbulent dynamics. The new equilibria and traveling waves that we present here form 
the backbone of this repertoire. Currently, a taxonomy of these myriad states eludes us, 
but emboldened by successes in applying periodic orbit theory to the simpler, warm- 
up Kuramoto-Sivashinsky problem (Christiansen et al 1997; Lan & Cvitanovic 2008; 
Cvitanovic et al. 2009), we are optimistic. Given a set of equilibria, the next step is to 
understand how the dynamics interconnects the neighborhoods of the invariant solutions 
discovered so far; a task that we address in Halcrow et al. (2009) which discusses their 
heteroclinic connections, and Gibson & Cvitanovic (2009) which discusses their periodic 
orbit solutions. 

The reader might rightfully wonder what the small-aspect periodic cells studied here 
have to do with physical plane Couette flow and wall-bounded shear flows in general, 
with large aspect ratios and physical spanwise-streamwise boundary conditions. Indeed, 
the outstanding issue that must be addressed in future work is the small-aspect cell 
periodicities imposed for computational efficiency. So far, most computations of invariant 
solutions have focused on spanwise-streamwise (axial-streamwise in case of the pipe flow) 
periodic cells barely large enough to allow for sustained turbulence. Such small cells 
introduce dynamical artifacts such as lack of structural stability and cell-size dependence 
of the sustained turbulence states. However, every solution that we find is also a solution 
of the infinite aspect-ratio problem, i.e., a solution whose finite [L X ,2,L Z ] cell tiles the 
infinite 3D plane Couette flow. As we saw in § 6, under a continuous variation of spanwise 
length L z such solutions come in continuous families whose fundamental wavelengths 
reflect the roll and streak instability scales observed in large-aspect systems such as 
figure 2. Here we can draw the inspiration from pattern-formation theory, where the most 
unstable wavelengths from a continuum of unstable solutions set the scales observed in 
simulations. 
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have independently determined the EQ 7 , EQ 8 equilibria.] 

REFERENCES 

Canuto, C, Hussaini, M. Y., Quarteroni, A. & Zang, T. A. 1988 Spectral Methods in 
Fluid Dynamics. Springer- Verlag. 

Cherhabili, A. & Ehrenstein, U. 1997 Finite-amplitude equilibrium states in plane Couette 
flow. J. Fluid Mech. 342, 159-177. 

Christiansen, F., Cvitanovic, P. k, Putkaradze, V. 1997 Spatio-temporal chaos in terms 
of unstable recurrent patterns. Nonlinearity 10, 55-70. 

Clever, R. M. & Busse, F. H. 1992 Three-dimensional convection in a horizontal layer sub- 
jected to constant shear. J. Fluid Mech. 234, 511-527. 



22 J. F. Gibson, J. Halcrow, and P. Cvitanovic 

Clever, R. M. k Busse, F. H. 1997 Tertiary and quaternary solutions for plane Couette flow. 

J. Fluid Mech. 344, 137-153. 
Cvitanovic, P., Davidchack, R. L. k Siminos, E. 2009 On state space geometry of the 

kuramoto-sivashinsky flow in a periodic domain. arXiv: 0709 . 2944, SIAM J. Appl. Dynam. 

Systems, to appear. 

Duguet, Y., Pringle, C. C. T. k Kerswell, R. R. 2008 Relative periodic orbits in transi- 
tional pipe flow. arXiv: 0807. 2580. 

Ehrenstein, U., Nag ATA, M. k Rincon, F. 2008 Two-dimensional nonlinear plane Poiseuille- 
Couette flow homotopy revisited. Phys. Fluids 20, 064103-1-4. 

Faisst, H. k Eckhardt, B. 2003 Traveling waves in pipe flow. Phys. Rev. Lett. 91, 224502. 

Frisch, U. 1996 Turbulence. Cambridge, UK: Cambridge University Press. 

Gibson, J. F. 2008a Channelflow: a spectral Navier-Stokes simulator in C++. Tech. Rep.. 
Georgia Inst, of Technology, Channelflow.org. 

Gibson, J. F. 20085 Movies of plane Couette. Tech. Rep.. Georgia Institute of Technology, 
ChaosBook . org/tutorials. 

Gibson, J. F. & Cvitanovic, P. 2009 Periodic orbits of plane Couette flow. In preparation. 

Gibson, J. F., Halcrow, J. k Cvitanovic, P. 2008 Visualizing the geometry of state-space 
in plane Couette flow. J. Fluid Mech. 611, 107-130, arXiv: 0705. 3957. 

Gilmore, R. k Letellier, C. 2007 The Symmetry of Chaos. Oxford: Oxford Univ. Press. 

Golubitsky, M. k Stewart, I. 2002 The symmetry perspective. Boston: Birkhauser. 

Halcrow, J. 2008 Geometry of turbulence: An exploration of the state-space of plane 
Couette flow. PhD thesis, School of Physics, Georgia Inst, of Technology, Atlanta, 
ChaosBook. org/projects/theses .html. 

Halcrow, J., Gibson, J. F., Cvitanovic, P. k Viswanath, D. 2009 Heteroclinic connections 
in plane Couette flow. J. Fluid Mech. 621, 365-376, arXiv: 0808. 1865. 

Hamilton, J. M., Kim, J. k Waleffe, F. 1995 Regeneration mechanisms of near-wall turbu- 
lence structures. J. Fluid Mech. 287, 317-348. 

Harter, W. G. 1993 Principles of Symmetry, Dynamics, and Spectroscopy. New York: Wiley. 

Hof, B., van Doorne, C. W. H., Westerweel, J., Nieuwstadt, F. T. M., Faisst, H., 
Eckhardt, B., Wedin, H., Kerswell, R. R. k Waleffe, F. 2004 Experimental obser- 
vation of nonlinear traveling waves in turbulent pipe flow. Science 305 (5690), 1594-1598, 
www . sciencemag . org/ cgi/reprint/305/5690/ 1594 . pdf . 

Hoyle, R. 2006 Pattern Formation: An Introduction to Methods. Cambridge: Cambridge Univ. 
Press. 

Itano, T. k Generalis, S. C. 2009 Hairpin vortex solution in planar Couette flow: A tapestry 

of knotted vortices. Phys. Rev. Lett. 102, 114501. 
Itano, T. k Toh, S. 2001 The dynamics of bursting process in wall turbulence. J. Phys. Soc. 

Japan 70, 701-714. 

J.E. Dennis, Jr., k Schnabel, R. B. 1996 Numerical Methods for Unconstrained Optimization 

and Nonlinear Equations. Philadelphia: SIAM. 
Jimenez, J., Kawahara, G., Simens, M. P., Nagata, M. k Shiba, M. 2005 Characterization 

of near- wall turbulence in terms of equilibrium and bursting solutions. Phys. Fluids 17, 

015105. 

Kim, H., Kline, S. k Reynolds, W. 1971 The production of turbulence near a smooth wall 

in a turbulent boundary layer. J. Fluid Mech. 50, 133-160. 
Kleiser, L. k Schumann, U. 1980 Treatment of incompressibility and boundary conditions 

in 3-D numerical spectral simulations of plane channel flows. In Proc. 3rd GAMM Conf. 

Numerical Methods in Fluid Mechanics (ed. E. Hirschel), pp. 165-173. GAMM, Viewweg, 

Braunschweig. 

LAN, Y. k Cvitanovic, P. 2008 Unstable recurrent patterns in Kuramoto-Sivashinsky dynam- 
ics. Phys. Rev. E 78, 026208, arXiv. org: 0804. 2474. 

Marsden, J. E. k Ratiu, T. S. 1999 Introduction to Mechanics and Symmetry. New York, 
NY: Springer- Verlag. 

Nagata, M. 1990 Three-dimensional finite- amplitude solutions in plane Couette flow: bifurca- 
tion from infinity. J. Fluid Mech. 217, 519-527. 

Nagata, M. 1997 Three-dimensional traveling- wave solutions in plane Couette flow. Phys. Rev. 
E 55, 2023-2025. 



Equilibria and traveling waves of plane Couette flow 



23 



Peyret, R. 2002 Spectral Methods for Incompressible Flows. Springer- Ver lag. 

Pringle, C. T. & Kerswell, R. R. 2007 Asymmetric, helical, and mirror-symmetric traveling 
waves in pipe flow. Phys. Rev. Lett. 99, 074502. 

Rincon, F. 2007 On the existence of two-dimensional nonlinear steady states in plane Couette 
flow. Phys. Fluids 19, 4105-+, arXiv : 0706 . 1165. 

Schmiegel, A. 1999 Transition to turbulence in linearly stable shear flows. PhD thesis, Philipps- 
Universitat Marburg, available on archiv.ub.uni-marburg.de/diss/z2000/0062. 

Schneider, T., Gibson, J., Lagha, M., Lillo, F. D. & Eckhardt, B. 2008 Laminar- 
turbulent boundary in plane Couette flow. Phys. Rev. E. 78, 037301, arXiv: 0805 . 1015. 

Skufca, J. D., Yorke, J. A. & Eckhardt, B. 2006 Edge of chaos in a parallel shear flow. 
Phys. Rev. Lett. 96 (17), 174101. 

Tuckerman, L. S. & Barkley, D. 2002 Symmetry breaking and turbulence in per- 
turbed plane Couette flow. Theoretical and Computational Fluid Dynamics 16, 91-97, 
arXiv : physics/0312051 . 

VlSWANATH, D. 2007 Recurrent motions within plane Couette turbulence. J. Fluid Mech. 580, 
339-358, arXiv: physics/0604062. 

VlSWANATH, D. 2008 The dynamics of transition to turbulence in plane Couette flow. In Math- 
ematics and Computation, a Contemporary View. The Abel Symposium 2006, Abel Sym- 
posia, vol. 3. Berlin: Springer- Verlag, arXiv : physics/0701337. 

Waleffe, F. 1990 Proposal for a self-sustaining mechanism in shear flows. Center for Turbulence 
Research, Stanford University/NASA Ames, unpublished preprint (1990). 

Waleffe, F. 1995 Hydrodynamic stability and turbulence: beyond transients to a self-sustaining 
process. Stud. Applied Math. 95, 319-343. 

Waleffe, F. 1997 On a self-sustaining process in shear flows. Phys. Fluids 9, 883-900. 

Waleffe, F. 1998 Three-dimensional coherent states in plane shear flows. Phys. Rev. Lett. 81, 
4140-4143. 

Waleffe, F. 2001 Exact coherent structures in channel flow. J. Fluid Mech. 435, 93-102. 
Waleffe, F. 2002 Exact coherent structures and their instabilities: Toward a dynamical- system 

theory of shear turbulence. In Proceedings of the International Symposium on "Dynamics 

and Statistics of Coherent Structures in Turbulence: Roles of Elementary Vortices" (ed. 

S. Kida), pp. 115-128. National Center of Sciences, Tokyo, Japan. 
Waleffe, F. 2003 Homotopy of exact coherent structures in plane shear flows. Phys. Fluids 

15, 1517-1543. 

Wang, J., Gibson, J. F. & Waleffe, F. 2007 Lower branch coherent states in shear flows: 

transition and control. Phys. Rev. Lett. 98 (20). 
Wedin, H. & Kerswell, R. R. 2004 Exact coherent structures in pipe flow: travelling wave 

solutions. J. Fluid Mech. 508, 333-371. 



